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Abstract. Since its commissioning in fall 1996, the 
stereoscopic system of Imaging Atmospheric Cherenkov 
Telescope (IACTs) of HEGRA with an energy threshold 
of 500 GeV, an angular resolution of 0.1° and an energy 
resolution of 20% per individual photon, and an energy 
flux sensitivity vY v at 1 TeV of 10~ n ergcuT 2 s" 1 (S/N 
= ho) for one hour of observation time has been used to 
monitor the BL Lac object Mkn 421 on a regular basis. 
In this letter, we report detailed temporal and spectral 
information about the TeV characteristics of Mkn 421 in 
1997 and 1998. We study the light curve, the shortest time 
scales of flux variability, the differential spectra on diurnal 
basis for several days with good 7-ray statistics and the 
time averaged energy spectrum. Special emphasis will be 
put on presenting the data taken during the world-wide 
April 1997 multiwavelength campaign. We compare the 
Mkn 421 results with the results obtained for the BL Lac 
object Mkn 501 and discuss possible implications for the 
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emission mechanism and the Diffuse Intergalactic Back- 
ground Radiation. 
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1. Introduction 

The advent of the Compton Gamma Ray Observatory and 
the EGRET instrument on board has led to the spectacu- 
lar discovery of intense high energy 7-ray radiation from as 
many as ~65 Active Galactic Nuclei (AGN) in the 0.1 to 
10 GeV energy range (Hartman et al. |1999| ). The recorded 
sources belong mainly to the blazar class, i.e. BL Lac ob- 
jects, Optically Violent Variables (OVVs), and highly po- 
larized quasars, all radio- loud sources which show variabil- 
ity on short time scales in most of the observed frequency 
bands. The power emitted in 7-radiation frequently dom- 
inates the power radiated by the source (see e.g. the spec- 
tral energy distributions shown by Fossati et al. 1998). 
It is widely believed, that the non-thermal 7-ray emis- 



sion is produced in a relativistic jet by a population of 
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electrons, emitting synchrotron radiation at longer wave- 
lengths and 7-rays in Inverse Compton (IC) processes of 
the highest energy electrons with lower energy seed pho- 



tons (see for recent reviews Coppi 1997; Sikora 1997; Ul- 
rich et al. 1997). The origin of the IC seed photons has 
not yet been established, debated possibilities include a 
target photon population dominated by lower energy syn- 
chrotron photons (the Synchrotron Self Compton mech- 



anism , see e.g. Bloom & Marsch er 1993; Ghisellini et al 



1996 ; Mastichiadis & Kirk 1997 ), or dominated by exter 



nal photons, e.g. by radiation from the nuclear continuum 
scattered or reproc essed in the broad-line region s (Sik ora, 
Begelmann & Rees 1994 ; Blandford & Levinson 1995|), o r 
by accretion disc photons (Dermer & Schlickeiser 1994 ). 
Although these possibilities have been studied in great 
detail, models where hadrons are the particles primarily 
accelerated (see e.g. Mannheim et al. 1993 ) can not be 
ruled out yet. Clearly, the study of the emission process 
could yield crucial clues regarding the structure of the jets 
leading ultimately to an understanding of the mechanisms 
of energy extraction operating in the surrounding of the 
central supermassive object. 

In this paper, we report on the 1997 and 1998 obser- 
vations of the BL Lac object Mkn 421 (redshift z= 0.031) 



with the IACT system of HEGRA (Daum et al. |1997|) . 
The system provides an energy threshold of 500 GeV, an 
angular resolution of 0.1° and an energy resolution of 20% 
for individual photons, and an energy flux sensitivity v F„ 
at 1 TeV of 10" 11 ergon^s" 1 (S/N = 5a) for one hour 
of observation time. It has been used to monitor Mkn 421 
on a regular basis ever since its commissioning in fall 1996 
to gain detailed temporal and spectral information. 

TeV observations yield information about the emit- 
ting particles with very high energies, complementary to 
the information from longer wavelengths. Due to the ex- 
treme conditions needed for the production of TeV pho- 
tons, TeV blazars could be extremely revealing labora- 
tories for the understanding of these sources in general. 
It should be pointed out, that the TeV radiation is ex- 
pected to be partially absorbed by the Diffuse Extra- 
galactic Backg round Radiation (DEBRA) ( e.g. G ould & 
Schreder 1965 ; Stecker, de Jager & Salomon 1992 ). While 
on the one hand the DEBRA extinction opens unique pos- 
sibilities to infer information about the DEBRA density 
in the largely unconstrained 0.5 to 100 /im wavelength re- 
gion (see Aharonian et al. 1999h| , called Paper 2 in the 
following), on the other hand it renders the understand- 
ing of the source more difficult (e.g. Coppi & Aharonian 
19991; Bednarek & Protheroe |l999j). 



The BL Lac Mkn 421 was the first extragalactic object 



to be discovered as a TeV emitter (Punch et al. 1992; Petry 
et al. |1996| ). The early HEGRA observations during the 
1994/1995 observation period showed, that the spectrum 
of Mkn 421 was si gnific antly softer than the Crab spec- 
trum (Petry et al. 1996 ). While the mean TeV flu x leve l 
of Mkn 421 is about 0.5 Crab units (Petry et al. |l996|) , 



flares exceeding five Crab units on timescales as short as 
15 minutes have been observed (Gaidos et al. 1996). 

The second BL Lac object to be discovered as TeV 
source (Quinn et al. |1996| ; Bradbury et al. 1997), Mkn 
501 (redshift z— 0.034), underwent a major outburst in 
X-ray and in t he TeV energy range during 1997 (see e.g. 
Protheroe et al. 1997[ ). Detailed accounts of the results ob- 
tained with the HEGRA IACTs are given in (Aharonian 
et al. |l999a , called Paper 1 in the following), Paper 2, and 
in (Aharonian et al. 1999c). The source emission showed 
strong variability with differential fluxes at 1 TeV from a 
fraction to 10 times the Crab flux. The HEGRA IACT sys- 
tem was used to determine diurnal differential spectra for 
more than 60 individual days. Surprisingly no correlation 
between absolute flux and spectral shape was found. The 
differential time averaged photon spectrum was well de- 
scribed by a power law with an exponential cutoff dN/dE 
cx (E/lTeV)" 1 ' 9 exp(-E/6.2TeV) up to energies of - 20 
TeV. Due to the similarity in redshifts of Mkn 421 and 
Mkn 501 (z = 0.031 and 0.034 respectively) the intergalac- 
tic extinction of the TeV photons caused by the DEBRA 
should be very similar for both objects and a comparison 
of both spectra is of utmost interest. 

The paper is organized as follows. In Sect. 0the IACT 
system of HEGRA is introduced, and in Sect.^] the data 
sample and the analysis techniques are described. The re- 
sults concerning the time averaged and the diurnal proper- 
ties, as well as the search for variability within individual 
nights are presented in Sect. [|. In Sect. [| the results are 
summarized and discussed. 



2. The HEGRA Cherenkov telescope system 

The HEGRA collaboration operates six imaging atmo- 
spheric Cherenkov telescopes located on the Roque de los 
Muchachos on the Canary island of La Palma, at 2200 m 
above sea level. A prototype telescope (CT1) started op- 
eration in 1992 and has undergone significant hardware 
upgrades since then (Mirzoyan et al. |1994 ; Rauterberg et 
al. 1995). This telescope continues to operate as an inde- 



pendent instrument and the Mkn 421 data taken with this 
telescope will be presented elsewhere. The stereoscopic 
system of Cherenkov telescopes consists of five telescopes 
(CT2 - CT6), and has been taking data since 1996, ini- 
tially with three and four telescopes, and since fall 1998 
as a complete five-telescope system. Four of the telescopes 
(CT2, CT4, CT5, CT6) are arranged in the corners of 
a square with roughly 100 m side length, and one tele- 
scope (CT3) is located in the center of the square. During 
1997 and 1998, when the data discussed in this paper were 
taken, CT2 was still used as stand alone detector. 

The telescopes have an 8.5 m 2 tessellated reflector, fo- 
cusing the Cherenkov light onto a camera with 271 Pho- 
tomultipliers (PMTs), covering a field of view of 4.3° in 
diameter. A telescope is triggered when the signal in at 
least two adjacent PMTs exceeds an amplitude of 10 (be- 
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fore June 1997) or 8 (after June 1997) photoelectrons; in 
order to trigger the IACT system and to initiate the read- 
out of data, at least two telescopes have to trigger simul- 
taneously. Typical trigger rates are in the 10-16 Hz range. 
The PMT signals are digitized and recorded by 120 MHz 
Flash- ADCs. The telescope hardware is described by Her- 



mann (1997); the trigger system and its performance are 
reviewed by Bulian et al. (1998); the calibration of the 



detector sensitivity is described in (FraB et al. 1997), and 



the determination of the correct telescope pointin g wit h 
35 arc sec accuracy is described in (Puhlhofer et al. 1997 ). 



3. Data sample and analysis tools 

The analysis described in the following is based on 165 
hours of data acquired between January 1st, 1997 and May 
27th, 1998 with Mkn 421 altitudes larger than 45° above 
the horizon. Only data taken with all four IACTs during 
excellent weather conditions, i.e. clear sky and humidity 
below 90%, have been used. Furthermore, data runs with 
a mean Cosmic Ray rate deviating by more than 21% from 
the expected value or with a mean width parameter per 
telescope deviating by more than 6% from the expected 
value were excluded from the analysis. 

The Monte Carlo simulations used for the data analy- 



sis ar e desc ribed in detail in (Hemberger 1998; Konopelko 
et al. |199SD . Changes of the trigger conditions, the mirror 
alignment and reflectivities as well as in the PM sensitiv- 
ities divide the data sample in eight phases of approxi- 
mately uniform detector performance. For each of these 
eight phases, the detector simulations have been adjusted 
to describe as closely as possible the hardware perfor- 
mance; hereby, information about the mirror alignment 
from point runs, information about the mirror reflectivi- 
ties and Photomultiplier sensitivities from the measured 
Cosmic Ray rates, and the measured single pixel trigger 
characteristics have been taken into account (see Paper 1, 
2 and references therein). The Monte Carlo fine tuning and 
the data quality have intensively been tested using data 
from Mkn 501, Mkn 421, and from the Crab nebula, al- 
lowing to effectively test all response functions relevant for 
the determination of 7-ray fluxes and 7-ray energy spectra 



(see Paper 1 & 2, and Hofmann 1997). The current status 
of the Crab analysis with the HEGRA telescope system is 



described in (Konopelko et al. 1998) 



The standard analysis methods described in detail in 
Paper 1 and 2 have been used. Thus, spectral studies have 
been performed with loose selection cuts, i.e. with a cut 
on the angular distance of the reconstructed direction 
from the Mkn 421 direction 6 2 < (0.22 ) 2 and a cut on 
the parameter "mean scaled width" with a cut value of 1.2. 
Loose cuts do not give an optimal signal to noise ratio but 
guarantee small systematic errors on the energy dependent 
cut efficiencies. The analysis uses an extended OFF-region 
for minimizing the statistical errors of the background es- 
timate. A ring segment (of 180° opening angle), from 0.3° 



to 0.7° distance from the camera center, at the opposite 
side of the ON-region, has been chosen. The energy re- 
construction required two IACTs within 200 m from the 
shower axis, each with more than 40 photoelectrons per 
image and a "distance" parameter of smaller than 1.7°. 
Additionally, a minimal stereo angle larger than 20° was 
required. 

The methods for reconstructing the 7-ray energy spec- 
tra, as well as the description of the main sources of sys- 
tematic uncertainties have been discussed in detail in Pa- 
per 1 and 2. The systematic errors can be divided into two 
parts, into an error on the absolute energy scale, and into 
additional errors on the shape of the spectrum (see Paper 
2). The error on the absolute energy scale amounts to 15% 
and is dominated by a 10% uncertainty of the detector 
sensitivity. Additional uncertainties arise mainly due to 
the incomplete knowledge of the atmospheric conditions 
during the observations. The error on the shape of the 
spectrum is considerable in the threshold energy region, 
i.e. at energies below 1 TeV. The error is dominated by 
uncertainties in the trigger behavior of the IACTs and by 
uncertainties of the reconstruction of 7-ray energy spec- 
tra based on Monte Carlo simulations at discrete zenith 
angles. Integral fluxes above a certain energy threshold 
are obtained by integrating the differential energy spectra 
above the threshold energy, rather than by simply scal- 
ing detection rates. By this means the integral fluxes can 
be computed without assuming a certain source energy 
spectrum. 

For the data up to zenith angles of 30°, comprising 
altogether 126 h of data, spectral results are presented 
for energies above 500 GeV; for the 39 h of data with 
zenith angles between 30° and 45°, an energy threshold of 
1 TeV has been chosen. If not stated otherwise, the results 
presented in the following are derived from the data with 
zenith angles below 30°. 

4. Experimental results 

4--1- Time averaged properties 

The total excess from Mkn 421 is shown in Fig. [I]. Using 
an angular cut of O 2 < (0.22 ) 2 the significance of the 
excess ((N on — N ff)/y/N on + Naff ) is 16er before im- 
age cuts and 26tr after the image cut. The total number 
of recorded 7-rays is about 8000 at hardware threshold. 
The time averaged energy spectrum, determined from the 
data with zenith angles smaller than 30° is shown in Fig. 
||. The error bars show the statistical errors, the hatched 
region shows the systematic error on the shape of the spec- 
trum described above. Additionally, there is a 15% error 
on the absolute energy scale. From the lowest energies of 
500 GeV to the highest energies detected at about 7 TeV, 
the spectrum can be described by a power law (Fig. |[ 
solid line) 

^ = (12.1 ±0.5 stat ±4.3 syst ) HT 12 x 
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Fig. 1. The Mkn 421 excess from 1997-1998. The full circles 
show the events in the ON-region, the dashed line the events 
in the OFF-region (background region) after a loose "mean 
scaled width" cut and a software threshold of at least 40 pho- 
toelectrons in two or more telescopes. 
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The 36% systematic error on the flux at 1 TeV takes into 
account the error on the shape of the spectrum and the 
15% error on the energy scale. The x 2 -value of the fit, 
computed with the statistical errors only, is 23.5 for 10 
degrees of freedom, corresponding to a chance probability 
for larger x 2 ~ vs ^ ues of 1%. Taking into account the sys- 
tematic errors on the shape of the spectrum, the result is 
surely compatible with a pure power law spectrum. Note 
however, that a power law with an exponential cutoff: 



dN/dE oc (E/TeV) 



(-2.5±0.4 Htat ) 



exp(-E/E ) 



(2) 



with En = 2. 



+2.0 
-0.9 



TeV, fits the data equally well 



(Fig. ||, dashed line). Here the x 2 - vaiue i a l s ° computed 
with the statistical errors only, is 16.8 for 9 degrees of 
freedom with a chance probability for larger x 2 -values of 
5%. The evaluation of the data from 30° to 45° zenith 
angles (39 h) gives similar results, although the statistics 
are very poor; the power law fit from 1 TeV to 5 TeV 
yields a differential spectral index of -3.16 ±0.3 s t a t- 

The statistics do not permit to decide whether the en- 
ergy spectrum changed from 1997 to 1998: for 1997, a 
differential spectral index of -3.28 ±0.20 s t a t is computed, 
and for 1998 it is -3.00 ± 0.05 8ta t- 

A possible correlation between absolute flux and spec- 
tral shape has been searched for by determining the time 
averaged energy spectrum for all days with a differential 
flux at 1 TeV, between 5- 10~ 12 cm^s^TeV" 1 
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Fig. 2. The time averaged 1997-1998 Mkn 421 energy spec- 
trum. Vertical errors bars indicate statistical errors. The 
hatched area gives the estimated systematic errors, except the 
15% uncertainty on the absolute energy scale. The solid line 
sfeows a power law fit, the dashed line a power law with an 
dsponjjntial girt off. Upper limits are at 2g.corrfj£jence^level. 
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Fig. 3. The figure shows the energy spectrum for all days 
with a time averaged differential flux at 1 TeV, between 
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Fig. 4. The 1997 (upper panel) and 1998 (lower panel) Mkn 
421 light curve, i.e., the integral fluxes above 1 TeV as a 
function of date. The upper light curve starts at MJD 50449 
(1.1.1997) and ends at MJD 50610 (11.6.1997). The lower light 
curve starts at MJD 50814 (1.1.1998) and ends at MJD 50975 
(11.6.1998). The full symbols show detections with at least 2a 
significance and the open symbols are detections with less than 
2a significance. The upper limits are 2a confidence level, (sta- 
tistical errors only; see text for systematic errors) 
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The first data sample contains 



and 15- 10 
30- 1CT 12 cm-V 
43 h of data, the second one contains 13 h of data. 
The two spectra are shown in Fig. || for all bins with 
reasonable statistical errors. A power law fit to the 
high state data gives the differential spectral index of 
-2.98± 0.09 s t a t, and for the low state data we obtain 
-3.02±0.12 sta t, not indicating any difference between the 
two spectra. 



4-2. Diurnal results 

Figure | shows the 1997 and 1998 Mkn 421 light curve. 
Only statistical errors are shown. The systematic error on 
the integral flux, dominated by the 15% uncertainty in the 
energy scale, amounts to a 40% relative uncertainty. Mkn 
421 was highly variable during 1997 and 1998; the fluxes 
varied from an emission of a fraction of one Crab level to 
more than one Crab level. Most of the flares seemed to last 
one or two days. In 1998, from April 17th to April 28th 
(MJD 50920 - MJD 50931) a high state lasting for approx- 
imately ten days was observed. This flare coincided with 
a world-wide multiwavelength campaign lasting ap proxi- 
matel y fro m April 10 to April 30 (Takahashi et al. |l99S; 



Urry |199SQ . A more detailed discussion of the HEGRA 
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Fig. 5. Two diurnal Mkn 421 energy spectra. Upper limits 
are at 2a confidence level, (statistical errors only; systematic 
errors as for Fig. m) 



data from this campaign is given in the next subsection. 
During the whole 1997-1998 period no flares with fluxes 
in excess of two Crab levels were observed. 

For days with long observation times or with strong 
source activity, diurnal differential spectra can be ob- 
tained. Two examples of spectra are shown in Fig. p§. The 
first spectrum was acquired during 1.6 hours and the sec- 
ond one during 3.4 hours. These diurnal spectra are statis- 
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Table 1. HEGRA results for Mkn 421 on diurnal basis (1997 and 1998). 
Only the statistical errors are given. The systematic error on the flux at 
1 TeV is about 40% and the systematic error on the spectral indices is 
about 0.10. 
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tics limited to energies below 4 TeV. In Table [j] the results 
of power law fits are summarized for all days with good 
statistics (3 sigma excess above 1 TeV). No highly sig- 
nificant evidence for spectral variability has been found; 
the diurnal differential spectral indices are statistically 
consistent with the mean index of -3.1. The two most 
significant deviations from the mean index were found 
for April 22nd/23rd, 1998 (MJD 50925/50926) and April 
27th/28th, 1998 (MJD 50930/50931), where the spectral 
indices of -2.71 and -2.82 deviated by 2.1cr and 1.8cr from 
the mean value, respectively. 

4.3. The data of April 1998 

As mentioned above, the data of April 1998 are of spe- 
cial interest since a world-wide multiwavelength campaign 
was performed with measurements in the radio, optical, 
and UV regimes, at X-Ray energies (RXTE, BEPPO Sax, 
ASCA), at GeV energies, and at TeV energies. 

In Figure |, the HEGRA integral fluxes above 500 GeV 
and above 1 TeV are shown (data up to 30° zenith angle 
in both plots). The flux varied by a factor of up to four. 
Note, that all days of the multiwavelength campaign are 
included in Table |l|, since the emission was high and the 
observation windows were large. 

In Fig. fj], the mean spectrum during the April 
campaign is shown. A power law fit yields a differential 
spectral index of -3.03 ±0.08 sta t ± 0.10 S y St . 
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Fig. 6. Integral HEGRA fluxes for Mkn 421 during the 
world-wide multiwavelength campaign for a threshold energy 
of 500 GeV in the upper panel and a threshold energy of 1 TeV 
in the lower panel. (The statistical errors of the 500 GeV data 
sample are larger than for the 1 TeV data sample, since events 
near the detector threshold of 500 GeV, which were taken 
usder larg er zenith angles ent e r the integral flux with large 
wgights.) 
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Fig. 7. The Mkn 421 spectrum during the April 1998 observa- 
tion period (17.-28.04.1998). The solid line shows a power law 
fit. Upper limits are at 2a confidence level, (statistical errors 
only; systematic errors as for Fig. 0) 
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The mean differential flux at 1 TeV was (23.5 ±l.l sta t 
±8.4 syst )10~ 12 cm -2 s -1 TeV -1 , which roughly equals the 
Crab flux at 1 TeV. The x 2 -value computed with the sta- 
tistical errors is 16.1 for 9 degrees of freedom, correspond- 
ing to a chance probability for larger x 2 -values of 6.5%. 



4- 4- Search for variability within individual nights 

Temporal variability within individual nights has been 
searched for using the integral fluxes above 500 GeV (data 
with zenith angles up to 30°) and above 1 TeV (data with 
zenith angles up to 45°). For the data of each night, the 
integral fluxes were determined with 15 min bins. For the 
first data selection the mean observation duration is 1.7 h 
per night, for the second one it is 25% higher, since data of 
larger zenith angles are included. Variability was searched 
for with a x 2 -analysis, fitting constant models to the in- 
tegral flux rates and determining the chance probabilities 
P c for the resulting x 2 -values. 

The results for the 1 TeV energy threshold data sam- 
ple are given in Fig. |[ On the left side the distribution 
of the -logio(P c )-values is shown. The chance probabilities 
distribute as expected in the absence of variability. On the 
side of small chance probabilities with -logio(-P c ) > 2, a 
marginally significant excess of three nights above an ex- 
pectation of 0.76 can be recognized. Taking into account 
that 76 nights have been studied, the chance probability 
for finding such three nights is computed to be 4%. The in- 
tegral fluxes of two nights are shown in Fig. || (right side). 
The lower panel shows one of the three nights with the 
most significant variability (May 21st/22nd, 1998, MJD 
50954/50955). The examples exemplify that due to lim- 
ited photon statistics, the search is only sensitive to strong 
flares. Also for the 500 GeV energy threshold data sample, 
no individual night with significant variability has been 
found. 

A preliminary analysis of Whipple data taken during 
the world-wide multiwavelength campaign indicated flux 
variability at energies above 2 TeV on a time scale of sev- 



result is consistent with earlier measurements of the differ- 
ential spectral index during periods of moderate flux, i.e. 
-3.6±1.0 sta t (at >1.2TeV) measured during the 1994/1995 
observation period with the HEGRA CT1 telescope (Petry 



eral hours (Maraschi et al. 1999). No simultaneous data 
was taken with the IACT system of HEGRA. 

5. Discussion 

In this paper we present detailed temporal and spectral 
information about Mkn 421 during 1997 and 1998. The 
mean 1997-1998 differential flux at 1 TeV was about one 
half of the flux of the Crab nebula. The light curve shows 
several distinct flares with mean durations of typically one 
or two days and fluxes at 1 TeV of approximately the Crab 
level. During the world-wide multiwavelength campaign in 
1998, we observed a 10 day phase of increased Mkn 421 
activity with a differential flux at 1 TeV of one Crab level. 

The time averaged energy spectrum follows to good 
approximation a power law with a differential spectral in- 
dex (500 GeV to 7 TeV) of -3.09±0.07 sta t ±0.10 syst . Our 



et al. 1996| ) and -2.92 ±0.22 sta t ±0.1 syst measured during 



the 1995/1996 observation period with the Whipple tele- 
scope (Zweerink et al. 1997). 

We did not find highly significant deviations from this 
mean spectrum, neither by dividing the data in a 1997 and 
a 1998 data sample, nor by dividing the data in a low and 
a high emission data set. Moreover, all diurnal spectra - 
some can be determined with an accuracy of 0.2 in the 
differential spectral index - are consistent with the mean 
spectrum. 

The Whipple group reported on an extraordinary Mkn 
421 flare observed on May 7th, 1996 with an integral TeV 
flux of more than 5 Crab units and a flux increase by 
a fac tor of ^2 within approximetaly 1 h (Gaidos et al. 
1996), which permitted to determine a differential spec- 



tral index during the flare with good statistical accuracy 
(Zweerink et al. 1997; Krennrich et al. 1999). An index 



of — 2.56 ±0.07 sta t ±0.1 syst indicated a harder spectrum 
in the flaring state than in the quiescent state. Whereas 
our Mkn 421 spectra are consistent with a spectral index, 
which is independent of time and absolute flux, just as in 
the case of Mkn 501 (Paper 1), the Whipple spectrum dur- 
ing the extraordinary flare appears to be harder. Never- 
theless, taking into account the statistical and systematic 
errors on the estimates of the photon power law index, the 
effect of spectral hardening during the very strong flare is 
only marginally significant. 

In our data we did not find strong evidence for flux- 
variability within individual nights, although the sensitiv- 
ity for this search was limited, due to short observation 
windows of approximately 2 h and the low levels of emis- 
sion. 

Due to the similar redshifts of Mkn 421 and Mkn 501, 
the comparison of their spectra promises insights into the 
intergalactic absorption by the DEBRA. In Fig. ||, the 
1997-1998 Mkn 421 spectrum is compared to the 1997 
Mkn 501 spectrum. During 1997 Mkn 501 was on average 
at 1 TeV more than 7 times brighter than Mkn 421 dur- 
ing 1997-1998. The Mkn 421 spectrum is much softer than 
the Mkn 501 spectrum. Between energies of one and sev- 
eral TeV the differential spectral index is -3.1 for Mkn 421 
and -2.2 for Mkn 501. Certainly, due to the similar red- 
shifts of both objects, the difference has to be attributed 
to different intrinsic spectra and not to differences in the 
extinction by the DEBRA. In the Mkn 501 spectrum we 
recorded an exponential cutoff at 6.2 TeV, which could 
well be the result of the intergalactic absorption by the 
DEBRA (see Paper 2 for a detailed discussion of possible 
origins of the exponential cutoff). As can be recognized 
from the fits to the Mkn 421 data of Eqs. [j] and |^, we can- 
not exclude an exponential cutoff for cutoff energies above 
~2.8 TeV; this question will further be studied in future 
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Fig. 8. The left side shows the distribution of the chance probabilities for finding for each night a stronger variability than 
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Fig. 9. The 1997 Mkn 501 and the 1997-1998 Mkn 421 spec- 
tral energy distributions. The hatched regions show the sys- 
tematic errors on the shape of the spectra (see text). Upper 
limits are at 2a confidence level. 



We anticipate that detailed modeling of the temporal 
and spectral characteristics of Mkn 421 during the 
multiwavelength campaign, combining the TeV data with 
informations from longer wavelengths, will permit to gain 
important constraints on the astrophysical conditions 
inside the jet. 
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